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We study the elastic properties of a single A/B copolymer chain with a specific sequence. We 
predict a rich structure in the force extension relations which can be addressed to the sequence. 
The variational method is introduced to probe local minima on the path of stretching and releasing. 
At given force, we find multiple configurations which are separated by energy barriers. A collapsed 
globular configuration consists of several domains which unravel cooperatively. Upon stretching, 
unfolding path shows stepwise pattern corresponding to the unfolding of each domain. While releas- 
^ ' ing, several cores can be created simultaneously in the middle of the chain resulting in a different 

path of collapse. 
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^ ■ I. INTRODUCTION 

o ■ 



Recent advances in nano manipulation allow direct access to single molecule elasticity using Atomic Force Mi- 
croscopy (AFM) 0-§]. AFM seems to be an excellent technique to study interactions in single chains by use of their 
elastic response. As a result of interplay between entropy, interactions between monomers and environments (such 
as solvents, other chains, external fields) the resulting force-extension relations allow deep insight into the nature of 
the interactions and conformations. In some polymer systems (including biopolymers and proteins), intrachain self 
Q : assembly produces secondary or tertiary structures and the elastic responses reflect this structural hierarchy 
^ ■ Furthermore, mechanically induced folding/unfolding leads a single molecule to undergo structural transitions foUow- 
, ing well defined specific paths without being frustrated. Especially in the case of proteins there is considerable hope 
^ ' to gain insight into folding mechanisms and design. 
^ ' In the present letter we study the conformation of random A/B copolymers and their responses to external forces. 
^-H As the response reveals information of the sequence and structure of the chain, this allows for "single chain force 
00 spectroscopy" . We use variational methods to provide some understanding of the elasticity of polymer with specific 
sequence. The present model corresponds also to the most simple "protein model" pO| ] when a polymer chain consists 
' of more than one constituents (here A and B monomers) which are arranged in a specific sequence along the backbone. 
The sequence carries relevant information since according to the arrangement of the monomers certain conformations 
are more preferred than others. 

Although the difficult nature of the problem suggests a detailed mathematical analysis (which is partly found 
■ below) , some physical scaling ideas help to understand the problem. Blob pictures and estimates of relevant length 
scales provide an instructive physical picture for homopolymer deformation in poor solvent condition. Here relevant 
I • length scale is the thermal blob size ^t- Each blob carries energy of k^T . A sharp structural transition occurs when 

the Pincus blob size ^p(= ksT/ f) reaches to this thermal blob size = ^t- 
C . Experimentally, two scenarios are available: either the displacement or the force are controlled [0-^. When the 
displacement is imposed and the force is measured, the unfolding of sequence of domains leads to repeating "saw 
tooth" pattern in force extension curve In constant force measurement, folding-unfolding transition appears at 

characteristic force, which is shown as "plateau" 1^. In the present study we minimize the free energy at given force 
(constant force measurement). As we will show later, we find multiple minima at a given force. These metastable 
states are created by the heterogeneity of the sequence. Some of these configurations are visited depending on the 
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' history of folding/unfolding. 



II. ANISOTROPIC VARIATIONAL METHOD 



Formally we consider a polymer chain consisting of N monomers of size b under tension which is described by 

H^H, + Hi + f-ifN-ri). (2.1) 

The first term corresponds to the elastic and connectivity properties of the polymer chain and is approximated by 
the Wiener measure [O]. The second term contains all interactions between monomers, especially, two and three 
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body interactions of the virial expansion are included here. The three body term is essential since depending on the 
solvent quality certain monomers may attract each other, which results in a (partial) collapse of the chain. The two 
body interactions can be attractive or repulsive depending on the type of the monomer pairs. The density of the 
collapsed globule is determined by balance between the three body repulsions and two body attractions. The last 
term is the elastic energy term which represents the external force field. The variational free energy is defined as 
usual by Fy = {H — Ht)o + Fq where (• • •)o stands for the average over the variational probability distribution: 

Pv{ri,..,fN) = Zy^cxp{-Ht{n,..,rN)}, (2.2) 

where Zy is the normalization constant satisfying J Py — 1 and Fq — —kBTlogZy. 

The presence of an external force breaks the spherical symmetry. We need, therefore, to distinguish the deformation 
in parallel and perpendicular direction with respect to the external force. In terms of two components in correlation 
function, we choose consequently the trial Hamiltonian Ht as follows 
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where G|| (j, /), G±{j, I) are correlation functions between j — I monomers in parallel and perpendicular direction with 
respect to the external force. By minimizing the trial free energy function with respect to G_l, G||, we obtain two sets 
of self-consistent equations. These self-consistent equations for both components are coupled, in contrast to systems 



with spherical symmetry as studied previously 12-14 



The square of the distance b{i,j) between two monomers i and j is related to G\\{i,j) by 

6|l = G|| (*, z) + G|| (j, j) - G|| - G|| (j, i). (2.4) 

Therefore, the square of the mean end to end distance (i?^(l,iV)) is determined by the quantity 6||(l,iV). The 
detailed derivation and solution of the self-consistent equation is shown elsewhere p^ . This procedure corresponds 
to the self-consistent one- loop approximation |16|. 



III. STRETCHING A HOMOPOLYMER IN A POOR SOLVENT 



As a test run we apply the variational calculation to calculate the elastic response of a homopolymer in a poor 
solvent. This problem has been studied theoretically (scaling ||l^, MC simulation [^), and experimentally p9[ . 

Using the present variational method, we confirmed the scaling picture of force-extension relations of homopolymer. 
The position of the plateau in each curve captures the characteristic force fc at each solvent condition and is in 
agreement with scaling picture. When we follow the local minima, we observe a hysteresis effect at the boundary of 
the first order phase transition. Around fc, we access to local minima with different extensions. This means that 
chain configuration becomes metastable. At the transition point, tadpole and ellipsoidal conformations coexist. 



IV. STRETCHING A HETEROPOLYMER IN A POOR SOLVENT 



Now we consider a random A/B polymer chain of length N. The interaction Hamiltonian for the pairwise interaction 
is given by 

N N 

i=i j=i 

where the two body interaction possesses a random nature, i.e., the interactions are chosen as a random N x N matrix 




Vij =vo + [a(fTj + aj) + x{o-i(Jj)]- (4.2) 

The sequence of monomers is described by variables (7^. = — 1 if monomer i is of type A (say, hydrophobic) and 
ai — 1 ii it is of type B (hydrophilic). To be more specific, a Flory like mixing parameter x — (vaa + wbb)/2 — vab 
is negative when similar monomers attract each other. The second virial coefficient vq < provides that average 
attraction between monomers. Obviously, the results will depend on the specific sequence of type A and type B 
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monomers. For technical reasons we assume that the average solvent quality is poor compared to the fluctuation due 
to randomness (|fo| < x)) therefore the thermodynamic ground state is always a globular state. Before a further 
discussion on the force extension properties for given sequences, we consider the physical effects of the randomness in 
our model to see what we can expect. 

V. STATISTICS 

To get a physical understanding for the problem we average over the disorder and assume a Gaussian distribution 
for the disorder part in the two body interaction, i.e., P(crj) = e"'^'/^'^. We obtain the first order correction in the 
effective binary second virial coefficient between arbitrary monomer pairs I — m: {v'^^}im = Vq — o? jx for ^ 7^ 

If the interactions between the similar monomers are equal (uaa = 'J^bb) then we have a = 0. We expect no effective 
change in the effective solvent quality. In our consideration, we assume that a ^ x.. The Flory mixing parameter x 
represents also the strength of disorder. Each thermal blob size within the collapse globular phase will be represented 
by the effective second virial coefficient Ueff instead of uq, which leads to = The average thermal 

blob size is ^™"<i°™ = ^homo^-|^ — a/|uo|). The mean square average of v'^^ is ^ a/^/N. For small disorder, for each 
sequence, the structural transition occurs at large scale. The variance in characteristic force /c reflects the size of 
disorder (~ a\fN). Although the mean net charge is zero, a typical random copolymer of size N has excess charge of 
order ±-\/iV. For some sequences, there might be an excess charge Y^di 4" 0- The characteristic force /c is changing 
accordingly due to the shift in an effective solvent quality. 

Depending on the mean square fluctuation, we divide the force-extension curve into three distinctive regime. For 
small extension (^p > / < k^T j ^^^^'^^^~) a collapsed globule deforms to ellipsoidal conformation. For large 

extension (^p < ^^^Uest^ j ^ fceT/^™^''''^*), all thermal blobs are linearly aligned. No conformational degree of 
freedom is left for rearrangement in thermal blobs. The blob size is then simply equal to Pincus blob size. In 
intermediate regime (^s™aUcst ^ ^ ^largest -j^ expect fluctuations in z{f) to depend on the system size and 
external force. 
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FIG. 1. Typical force-extension curve which show multiple steps for four different sequences with x/l'f'ol = 0.5. The hysteresis 
during stretch-release cycle is shown at each curve. The sequence of each polymer is 

a) AAAAB BBBAB BBBBB AABBA BABAA BAABA BAAAB BABAA 

b) BABAA BABBB BBAAA ABBAB BAAAB ABBAA ABBAB BBAAA 

c) ABAAB BAAAA AAAAB BAABB BBAAB BBBAB ABBBB BABAA 

d) AABAB AABBB AABAB BABBB AAABB BABAB ABBAA BABBA 
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VI. FORMATION OF DOMAINS: BLOCKS 



Due to the heterogeneity of the sequences, in average, the structural transitions are smoothened. However, unfolding 
of individual chains shows unique elastic response reflecting the information in the interactions along the sequence. 
Within a collapsed globular phase, a local sequence of "n" monomers has an excess charge with typical fluctuations of 
size \/n. This creates more hydrophobic (type A dominant) cores. Their sizes are determined by the disorder (in local 
sequence) and their positions along the chain are random. In fig. we show typical examples for the force-extension 
relationships which show multiple steps for four different types of sequences with x/|wo| = 0.5. 

Multiple plateaus indicate that the force induced unfolding from the globular to open-string structure occur via 
intermediate structural forms. Each plateau in force-extension relation corresponds to the unfolding of a part of 
the chain cooperatively which defines a "domain" . The structural changes are achieved by abrupt unraveling of 
the corresponding domain. It is clearly visible that the curves depend on the sequence. For example, curve d) in 
fig. |l| corresponds to the most "alternating" sequence (alternating copolymer). Therefore, the overall poor solvent 
dominates, and the conformation is a single globule. In the situations of curves a) and d), several A and B blocks 
are separated by "non-blocky" units. These allow for a more structured phase space, the force-extension relation is 
richer. Moreover sequence c) has a net charge 2 (zero charge for a,b and d). 




monomer index / monomer index / 

FIG. 2. The inverse correlation function h\\ (i + 5, i) during the stretch-release cycle. The chain conformation choose different 
paths for folding and unfolding procedure. The correlation 6||(i,j) of some part of chain increase/decrease cooperatively. 
Different symbols indicate applied force in units of ksT . 

In fig. H the correlation function of a specific sequence (c in fig. |l|) is shown during the stretch-release cycle. The 
regions of small values of 6||(i, j) represents the globular phase. As the force increases, a part of the chain is pulled 
out from one end leading to a tadpole like conformation. Upon further increase of the force, the remainders in 
the globular phases are reduced. Then, the unfolding of the globule occurs in stepwise patterns, which have also 
been predicted for polyelectrolyte necklace chains p3| . Different configurations are visited while the chain releases. 
The more hydrophobic parts assemble together at the beginning of the release, the necklace like conformations are 
formed. These clusters are growing and finally merge to a single globule. The fig. ||also shows that upon releasing the 
structural change from necklace to globule occurs at smaller characteristic force than upon stretching. The hysteresis 
on force extension curves is determined by the energy barriers between two coexisting configurations following a force 
induced path. The time for relaxation is exponentially increasing with respect to the energy barrier between 
the two states, ty. ^ ex.Y'{AE^ / k-oT) . These local minima created by quenched disorder is stable with respect to the 
fiuctuation . The hysteresis appears in the experimental situation with finite pulling speed . It has been already 
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seen that for pulling speeds faster than tr , the rearrangement of monomers at given displacement is not followed by 
unfolding of each domain. However, a complete analysis of the dynamical aspects and the influence of the pulling 
velocity will be given in a separate paper fl^ . 





FIG. 3. Blob interpretation of stretching a random copolymer 

The stretching without relaxation is illustrated in flg. ^ The external force propagates from ends to the center of 
the chain. The extension is controlled by smaller thermal blob size along the path of force propagation. Variations 
in thermal blob size due to the disorder along the chain are expected. Each blob carries the thermal energy fceT 
and the variation in size of blob comes from the variation in local sequence. Applying external force unravels each 
thermal blob one after another when the force reaches the order of = (i=A,B,C,...). The force needed to pull 
out blob B is large enough to pull blobs C and D (since £,c,^d < £,b)- Therefore, as soon as blob B is released, there 
is no further increase of the force required to pull C and D. This means that the B-C-D blobs respond to the external 
force cooperatively. In similar way, E-F-G form another domain. After being pulled out, the size of each blob should 
be rescaled to the size of Pincus blob at given force. 

In thermodynamic equilibrium at given force /o, more than one configuration with different extension can be 
accessed with some probability. (For example, 1 and 2 in fig. |] d). The ratio between the average life time of two 
configurations is P{Ei)/P{E2)- This is related to the free energy difference AF = ~kBTln{P{Ei)/P{E2))- From 
the force extension curve, the free energy difference can be easily obtained: AF = Fi — F2 — fodzi-2- Estimation of 

energy barrier along the path l-l'-2'-2 (shown in Fig. |^ d) is AEl_2 = f{z)dz. 



VII. CONCLUSION 



We investigated the elastic properties of a single chain with arbitrary sequences using variational method. The 
essential property governing the force-extension curve is domain-wise unfolding from globular to open string con- 
formation. The characteristic force related to each domain is captured as a "plateau" in the force-extension curve. 
In this variational approach, other types of interactions (including long range interactions) and other properties of 
the chain (e.g., stiffness) can be easyly incorporated. When long range interactions are dominant, domains becomes 
correlated. This effect will smoothen out the sharp transition due to the unfolding of each independent domain. We 
discuss the role of long range interactions separately . Proteins could be also treated in similar way by assigning 
special sequence using binary interaction matrix and by adding more types of monomers which may correspond to 
the amino acids. The accuracy of the method is tested in comparison with Monte Carlo simulation ||25[|. 
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